Life time of superconductive state in long overlap Josephson junction 
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The computer simulations of fluctuational dynamics of the long overlap Josephson junction in 
the frame of the sine-Gordon model with a white noise source have been performed. It has been 
demonstrated that for the case of constant critical current density the mean life time (MLT) of 
superconductive state increases with increasing the junction's length and for homogeneous bias 
current distribution MLT tends to a constant, while for inhomogeneous current distribution MLT 
quickly decreases after approaching of a few Josephson lengths. The mean voltage versus junction 
length behaves inversely in comparison with MLT. 



At the present time in the area of quantum calcula- 
tions serious hopes are connected with the possibility to 
create qubits on the basis of point 1], and distributed 
Josephson junctions Q The main advantage of qubits 
based on Josephson junctions is the relative simplicity of 
manufacturing and integrating in one circuit in compar- 
ison with other possible realizations of qubits 3J. It is 
necessary to point, that at the present time all Joseph- 
son junctions are manufactured with the use of electron- 
beam lithography 3,Q and can always be considered as 
distributed. One of the most fundamental problems of 
all quantum calculations is connected with the destruc- 
tion of entangled states of qubits through the interaction 
with the environment. This event is characterized by 
the time of decoherence, the maximal increase of which 
is needed for the design of quantum computers. At the 
same time, due to mathematical difficulties, even in the 
case of a long Josephson junction in the presence of ther- 
mal fluctuations, the life time of superconductive state 
is not studied well enough (see, e.g., 0-113). Mostly, 
annular and infinitely long junctions have been consid- 
ered, but it has been understood, that for short junc- 
tions the escape of the whole string over the potential 
barrier occurs, while for long junctions the creation of 
kink-antikink pairs is the main mechanism of the escape 
process. In our recent works, dedicated to investigation 
of noise-induced errors in Josephson junctions 
only the model of point junction was considered, so no 
dependences between fluctuational properties and junc- 
tion's length were studied. Therefore, at the present time 
an important problem of investigation of superconductive 
state life time with the aim to improve noise immunity of 
long Josephson junctions is not yet solved. In this Letter 
we present results of numerical investigation of supercon- 
ductive state life time based on the computer simulation 
of the sine-Gordon equation with white noise source. 

In the frame of the resistive McCumber-Stewart model 
[l3| the phase difference of the order parameter <p(x,t) 
of long Josephson junction in the overlap geometry (see 
Fig. 1) is described by the sine-Gordon equation: 

with the following boundary conditions: 



i(x)-sm((p) + i f (x,t), (1) 
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FIG. 1: The structure of distributed Josephson junction of 
the "overlap" geometry. 
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Here the time and the space are normalized to the inverse 
characteristic frequency uj~ 1 and the Josephson length 
Aj, respectively, (3 = 1/a 2 is the McCumber-Stewart 
parameter, a is the damping, i(x) is the bias current 
density, normalized to the critical current density of the 
junction, if(x,t) is the fluctuational current density, T 
is the normalized magnetic field, L — l/Xj is the di- 
mensionless length of the junction. In the case where 
the fluctuations are treated as white Gaussian noise with 
zero mean, and the critical current density is fixed, its 
correlation function is: 



where 



(i f (x,t)if(x',t')) = 2 1 8{x-x')5{t-t'), (3) 



7 = It/{Jc\j) (4) 



is the dimensionless noise intensity, J c is the critical cur- 
rent density of the junction, It = 2ekT/h is the thermal 
current, e is the electron charge, Ti is the Planck constant, 
k is the Boltzmann constant and T is the temperature. 
It should be noted, that in [g the total critical current 
I c = f J c (x)dx (here I — LXj) was supposed to be fixed, 
and the noise intensity in that case was proportional to 



the dimcnsionlcss length of the junction L. In the frame 
of the present paper we consider more interesting (from 
the practical point of view) case where the critical current 
density is held constant and the critical current increases 
with the junction's length. 

Numerical solution of the sine- Gordon equation JTJ 
with boundary conditions J3J) has been performed on the 
basis of implicit finite-difference scheme |l 11 with the ac- 
count of the white noise source. Typical values of dis- 
cretization steps are Ax = At = 0.1 — 0.02, number of 
realizations are R = 1000 - 10000. 

In the frame of the present paper we will consider and 
compare two cases of bias current distribution: homoge- 
neous i(x) — in and inhomogcneous one, characteristic 
for a superconductive thin film 13], 15]: 



i(x) 



i L 



~K\Jx(L — x) 



(5) 



For the junctions with small capacitance (3 <C 1 (large 
damping a 3> 1), the mean life time (MLT) of supercon- 
ductive state is, physically, the mean time till the gener- 
ation of noise-induced voltage pulse. The MLT is defined 
as the mean time of phase if existence in the considered 
interval 

[13, El, where W(ip,x,t) is the probability den- 
sity, and P(t) is the probability that the phase is located 
in the initial potential well: 



t = I tw{t)dt = 
'o Jo 



j(t) = 



dP(t) 
dt : 



P(t)dt, 



W(tp, x, t)dipdx. 



(6) 



First, let us consider the case of homogeneous bias 
current distribution i(x) — in. The MLT of supercon- 
ductive state versus length of a Josephson junction, cal- 
culated numerically on the basis of definition JHJ with 
the delta-shaped initial distribution at ipo = arcsin(io), 
is presented in Fig. 2. As it is seen, the MLT of super- 
conductive state for the case of constant critical current 
density rises for small lengths L < 1 and reaches a con- 
stant for L > 5. Such behavior is explained by the ex- 
istence of two different mechanisms of the string escape 
over the potential barrier || (see also 0-0] ) : it may es- 
cape either as a whole, or by forming kink-antikink pairs. 
For L < 1 with the increase of the string length, by the 
interaction of elementary parts of the string with their 
neighbors, the phase is stronger kept in the initial poten- 
tial well which leads to the increase of the life time. Even 
if the kink-antikink pair is fluctuationally formed, but the 
string is too short, the attraction between the kink and 
antikink outweighs the driving force i and the incipient 
nucleus collapses returning the string to the initial po- 
tential well. If, however, the string is long enough, then 
the driving force predominates and the nucleus expands, 
speeding up the escape process. These mechanisms are 



illustrated in Fig. [21 while for L < 1 the first mecha- 
nism is predominated, for L > 5 the second one leads to 
the saturation of MLT. It is important to mention, that 
for the considered overlap junction (see boundary con- 
ditions (J2J) not only the creation of kink-antikink pairs 
(curve 4), but also creation of one (curve 2), two (curve 
3), or even more kinks is possible. As it is seen from Fig 
13 the constant of saturation depends on thermal noise 
intensity 7 (decreases with the increase of 7), and bias 
current density i (also decreases with the increase of i, 
because of decrease of the potential barrier height). 
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FIG. 2: The MLT of superconductive state of long Josephson 
junction for the case of homogeneous bias current distribu- 
tion with r = 0, P = 0.01: squares - in = 0.5, 7 = 0.3; circles 
- in — 0.7, 7 = 0.3; triangles - in — 0.9, 7 = 0.3; solid line - 
7 = 0.2, i = 0.7; dashed line - 7 = 0.4, i = 0.7. 
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FIG. 3: Phase evolution in the long Josephson junction de- 
pending on its length: 1 - L = 2, 2 - L = 5, 3 - L = 20, 
t = ti, 4 - L = 20, t = t 2 > h. 

It is interesting to check the limiting transition of a 
long junction to a point one for small junction lengths 
L C 1. For a point junction in the case of large damping 
/3 C 1 the MLT of superconductive state was found an- 
alytically 0] for arbitrary value of noise intensity. That 
formula, however, was obtained for the case of constant 
critical current I c , and the noise intensity in that case 



2 




0-1 , , , , , , , ,— 1 , , , , , , , r- 

0.0 0.5 1.0 1.5 2.0 5 10 15 20 

L, arb. un. L, arb. un. 



FIG. 4: Comparison between numerical simulations of MLT 
(symbols) and formula £J (lines) for the case of homo- 
geneous bias current distribution with V — 0, [3 — 0.01, 
io —0.7; circles and solid line - 7 = 1, squares and dashed 
line - 7 = 0.3, triangles and dotted line - 7 = 0.2. 

was: 7 S = It/Ic, where I c — J c XjL. Substituting 7 S 
as j s — j/L (where 7 is given by J3J), one can get the 
following formula, that should describe MLT of supercon- 
ductive state for the case of a short Josephson junction 
with constant critical current density: 

T = - J f ' e -(cosx+tx)L/ 1 I e (™sV+*>f)Lh dlpdx+ 
p(p2 poo \ 

+ e {cos * +l * )Lh dv I e- {cos * + ^ )Lh dtp Y (7) 
Jip 1 J tp 2 ) 

where ipo - coordinate of initial delta-shaped distribu- 
tion, 1^1,2 - boundaries of the interval, which define the 
initial potential well. Indeed, as it is seen from Fig. 4, 
the results of numerical simulations of a long junction 
perfectly coincide with the formula J7J) not only in the 
limiting case L -C 1, but even up to L ~ 1 and for L > 1 
the coincidence is better for larger noise intensity. 

In real situations it is hard to achieve the homogeneous 
bias current distribution. Since a long overlap Josephson 
junction is made on the basis of thin superconductive 
films, the bias current distribution may be governed by 
the formula © 1713.11"") for the case when the junction 
width is much smaller than its length. Such bias current 
profile significantly changes the picture of the life time of 
superconductive state dependence on junctions's length. 
Results of study of the MLT in the case of inhomogeneous 
bias current density are presented in Fig. [SJ 

The main feature of such dependences is the existence 
of strong maximum for the lengths L « 5. This effect can 
also be explained by the existence of two mechanisms of 
the string escape. It is obvious, that mechanism of MLT 
increase for small lengths is absolutely the same as for 
the case of constant bias current distribution (see Fig. 
EJ. The main difference between the cases of homoge- 
neous and inhomogeneous bias current distribution is the 
fact, that in inhomogeneous case there are unstable areas 



FIG. 5: The MLT of superconductive state of long Joseph- 
son junction for the case of inhomogeneous critical current 
density with r = 0, ft = 0.01, 7 = 0.4 and: circles - io — 0.5, 
squares - io = 0.6, triangles - io = 0.7. 

near the edges of a junction (x w 0, L) because there the 
bias current is larger than the critical current according 
to the formula ©. These areas work as generators of 
" kink-antikink" pairs. Due to the unstable character of 
a potential profile near the edges, amount of those pairs 
is significantly larger than in the junction with homoge- 
neous bias current distribution. So, the influence of the 
second mechanism on the life time for the inhomogeneous 
bias current distribution is larger in comparison with the 
homogeneous one. But the lengths of "switching on" of 
this mechanism are the same as for the homogeneous case 
(see Fig. 01 . 
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FIG. 6: Comparison of MLT of superconductive state of long 
Josephson junction for different cases of bias current distri- 
bution with r = 0, f3 = 0.01, 7 = 0.3, i = 0.7: circles - 
inhomogeneous distribution, squares - homogeneous distri- 
bution. 

On one hand, to experimentally observe the predicted 
above effect of non-monotonous dependence of the MLT 
of superconductive state versus junction's length, one 
needs to perform advanced measurements. On the other 
hand, it is now believed that the use of more ho- 
mogeneous bias current distributions may increase the 
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power and decrease the lincwidth of practical Flux-Flow 
Oscillators (FFOs), based on long Josephson junctions. 
There are attempts to increase the homogeneity by the 
increase either the junction's width or the width of the 
idle regions (see [lq for details) and the check of the ho- 
mogeneity of the current profile is only performed by the 
measurements of the linewidth. This takes a long time 
and there are many other different factors which affect 
the linewidth. So, a reliable tool to detect a homogeneity 
of bias current distribution in a particular FFO design is 
needed. To this end, the measurement of a voltage-length 
characteristic for a fixed value of bias current (slightly 
smaller than the critical current) can be recommended 
as such a tool. 
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FIG. 7: Comparison of mean voltage V(L) for the differ- 
ent cases of bias current distribution with io = 0.9, T = 0, 
(5 = 0.01; 7 = 0.05: circles - inhomogeneous case, squares - 
homogeneous case; 7 = 0.1: solid line - inhomogeneous case, 
dashed line - homogeneous case. 

The mean noise-induced voltage versus junction's 
length demonstrates the "inverse" behavior in compar- 
ison with the MLT (see Fig. |7J). For the case of constant 
current density and homogeneous bias current distribu- 
tion the mean voltage decreases for small lengths L < 1 
and reaches a constant for L > 5. For the case of in- 
homogeneous current distribution the minimum of the 
mean voltage versus length is observed. As one can see 
from Fig. |7| the voltages for homogeneous and inhomo- 
geneous cases differ exponentially in the small noise limit 
for large lengths of Josephson junctions. If, with increase 
of the junction's width, the bias current distribution will 
be more and more close to the homogeneous one, it can 
be clearly seen in V(L) dependence. 

In this work numerical and qualitative analysis of fluc- 
tuational dynamics of long Josephson junctions is pre- 
sented. It has been demonstrated that for the case of con- 
stant critical current density the mean life time (MLT) of 
superconductive state increases with increasing the junc- 
tions 's length and for homogeneous bias current distri- 
bution MLT tends to a constant for L > 5. However, 
for inhomogeneous current distribution MLT quickly de- 



creases after approaching a maximum for lengths around 
L ~ 5. Therefore, from the fluctuational stability point 
of view, there is no reason to increase the length of a long 
Josephson junction more than L « 5, excepting the cases, 
where increasing of the junction length can improve other 
useful properties of Josephson electronic devices. It has 
also been demonstrated, that the measure of the homo- 
geneity of bias current distribution may be performed on 
the basis of voltage versus length characteristics for bias 
current, smaller than the critical one: for junction lengths 
larger than 5Aj, the voltages for bias current distribution 
(|5|) and the homogeneous one differ exponentially in the 
limit of small noise intensity. 
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